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Abstract. We discuss the ground state magnetic phase diagram of the Hubbard model off half filling 
within the dynamical mean-field theory. The effective single-impurity Anderson model is solved by Wil- 
son's numerical renormalization group calculations, adapted to symmetry broken phases. We find a phase 
separated, antiferromagnetic state up to a critical doping for small and intermediate values of U, but could 
not stabilise a Neel state for large U and finite doping. At very large U, the phase diagram exhibits an 
island with a ferromagnetic ground state. Spectral properties in the ordered phases are discussed. 
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1 Introduction 

Originally proposed for the description of ferromagnetism 
in transition metals, the Hubbard model [1] 



H 



^E^ 



is the simplest model to describe the interplay between 
delocalization or band formation in solids and the local 
Coulomb correlations. Despite its simplicity, the phase di- 
agram of the Hubbard model (1) reveals a surprising rich- 
ness. One finds Mott-Hubbard type metal-insulator tran- 
sitions [2], antiferromagnetism [3], ferromagnetism [4] and 
incommensurate magnetic phases [5]. More recently, the 
Hubbard model has also become one of the most promis- 
ing candidates to describe the low-energy properties and 
possibly the superconductivity in the high- Tc cuprates [6] . 

The Hubbard model at half filling (n) — 1 has been in- 
vestigated thoroughly using various approximate and ex- 
act techniques [7,8] and its properties are understood to a 
large extent. Off half filling, the model is well understood 
for d = 1 but the situation is less clear in dimensions o? > 1 . 
Basically the only rigorous result is due to Nagaoka [9], 
who proved that a ferromagnetic ground state is possible 
under certain conditions. 

The introduction of the limit D oo [10] in princi- 
ple allows to solve models like the Hubbard model exactly 
without loosing the competition between kinetic energy 
and local Coulomb repulsion [11]. This surprising insight 
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has subsequently triggered a large amount of investiga- 
tions of the infinite dimensional Hubbard model [11]. In 
addition, D ~^ oo turned out to be a reasonable start- 
ing point for weak-coupling expansions [13,14]. Within 
this approach, in addition to the expected magnetic or- 
der for a bipartite lattice, phase separation was found for 
the whole region of the magnetically ordered phase [14]. 
Since this result is based on a weak-coupling expansion, 
it is far from clear whether it holds for finite values of 
the interaction U as well. Results from a numerically ex- 
act solution of the Hubbard model in Z? = oo based on 
Quantum Monte-Carlo simulations for example showed no 
evidence for phase separation [15], but these calculations 
were done in the paramagnetic phase and at finite, com- 
paratively high temperatures. 

The question whether phase separation in the Hubbard 
model occurs in a certain parameter regime is of some im- 
portance for two reasons. First, from a model theoretical 
point of view, it is of course interesting to explore the sta- 
bility of the different possible ordered phases which might 
be unstable with respect to phase separation. Second, a 
vicinity to phase separation has been discussed as one of 
the possible ingredients to the superconductivity in the 
high- Tc cuprates [16,17]. Moreover, a tendency towards 
phase separation together with the long-range part of the 
Coulomb interaction may in principle lead to charge or- 
dered states such as stripe-phases. 

Phase separation has long been predicted [16] and in- 
deed been observed for the t-J model in Z? = 1, 2 [18]. 
Since the t-J model for vanishing J is connected to the 
Hubbard model in the limit U/t oo [19], additional in- 
formation about phase separation in the strong coupling 
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limit could thus be obtained. The critical J, below which 
phase separation vanishes in the t-J model is typically 
of the order of t [18]; this seems to rule out phase sep- 
aration in the Hubbard model at least for C//t ^ oo in 
D = 2. However, the situation has not been clarified yet, 
because direct inspection of the Hubbard model in the 
limit U/t^oo leads to contradictory results [20,21]. 

The results for the 2D Hubbard model for finite U 
available so far have not revealed signs for phase separa- 
tion [18,22]. However, these results are typically based on 
Quantum Monte-Carlo or related techniques, which have 
severe problems in the interesting parameter regime close 
to half filling and at very low temperatures. Consequently, 
one either has to restrict oneself to rather small system 
sizes [18] or use further approximations [22]. Thus, to our 
present knowledge a detailed study of the ground state 
phase diagram of the Hubbard model in the thermody- 
namic limit and in the vicinity of half filling, compris- 
ing weak, intermediate and strong coupling within a non- 
perturbative approach is not available. 

Such an approach is provided by the limit D oo, 
which allows for in principle exact calculations in the ther- 
modynamic limit for all model parameters, even at T = 0. 
The price one pays is the neglect of non-local dynamics, 
which of course is most severe for D < 2. Nevertheless, 
the theory can give valuable information about whether 
phase separation is possible at all. 

The paper is organized as follows. In section 2, we give 
a brief description of the solution of the Hubbard model in 
the limit D ^ oo. Results for the phase diagram and the 
dynamics in the different phases are presented in section 
3. The main conclusions of the paper are summarized in 
section 4. 



2 Theoretical background 

2.1 General remarks 

The dynamical mean-field theory (DMFT) to exactly solve 
the Hubbard model in the limit D — > oo is based on 
the work by Metzner and VoUhardt [10] and is by now 
well-established [11]. The basic ingredient is that for D 
oo the proper single-particle self energy z) becomes 

purely local or momentum independent, i.e. Z'(k, z) ^-^^ 
S{z) [10, 12]. This can be used to map the Hubbard model 
(1) onto an equivalent quantum impurity problem supple- 
mented by a self-consistency condition [11]. The remaining 
problem (the solution of a quantum impurity model) is, 
however, highly nontrivial. Several approximate and nu- 
merically exact techniques are currently available [11,23]. 

Most of these methods cannot access T — > or are 
restricted to the weak-coupling regime of the Hubbard 
model. The most reliable technique to solve the quantum 
impurity problem for all interaction strenghts U and fill- 
ings n at T = and low T is the numerical renormalization 
group (NRG) [23,24]. Originally, this method was set up 
to treat the paramagnetic problem only [24], but recent 
extensions have shown that calculations with a symmetry 



breaking field are possible with a similar level of accu- 
racy, too [25,26]. Hence we are able to study magnetically 
ordered phases directly at T = 0. 

In contrast to the standard NRG, a more refined ap- 
proach has to be used to calculate dynamical quantities 
in the presence of a magnetic field. This has first been 
noted by Hofstetter, who observed discrepancies in the 
magnetization calculated from the spectral functions and 
the ground state occupation numbers [26] . To resolve this 
problem, he proposed a modification of the standard method 
[27] to calculate the spectral function. A more detailed dis- 
cussion of this technical point and its physical background 
will be presented elsewhere. 

There are in principle two ways to determine the phase 
boundary between the paramagnetic and a magnetically 
ordered state. First, one can calculate the susceptibility 
corresponding to the anticipated order and look for a di- 
vergence. Second, one can allow for a proper symmetry 
breaking in the one-particle Green function and search 
for the region in parameter space where a solution with 
broken symmetry becomes stable. Especially for T = 
the first method is rather cumbersome in general and, by 
construction, also makes no statement about the thermo- 
dynamic stability of phases beyond the critical point. 

We thus use the second approach as our method of 
choice. However, this prohibits the search for incommen- 
surate phases, because only broken symmetries with a 
commensurate wave vector can be implemented that way. 
Since we arc interested mainly in standard Neel type anti- 
ferromagnetic order, the proper way is to introduce an AB- 
lattice structure and allow for difi^erent sublattice magne- 
tizations. The resulting Green function then becomes a 
2x2 matrix, which within the DMFT has the form [11] 

where k is a vector in the magnetic Brillouin zone (MBZ). 
For the Neel state on an AB-lattice a further simplifica- 
tion arises from the symmetry S^{z) = Sf{z). For the 
calculation this means that we do not have to solve inde- 
pendent quantum impurity models for the two sublattices, 
but only one for say sublattice A. 

But how to control the filling if the homogeneous so- 
lution (homogeneous concerning the charge distribution) 
turns out to be unstable towards phase separation? Fixing 
the chemical potential fj, is not sufficient as the system will 
be driven to a filling corresponding to a stable solution, 
such as n = 1. To enforce a metastablc state with finite 
doping we adopt a procedure which has been already used 
in calculations for the half-filled Hubbard model in a ho- 
mogenous magnetic field [28]. The schematic fiow diagram 
of the resulting DMFT self-consistency cycle is shown in 
Fig. 1. Starting from a paramagnetic solution for the de- 
sired doping, a homogenous or staggered symmetry break- 
ing is introduced and the corresponding effective medium, 
Fo-(cj), for the DMFT cycle is determined [11]. Keeping 
the medium r^{w) fixed, one now varies the on-site en- 
ergy of the effective SIAM until the desired doping has 
been reached. This result is used to obtain a new efi^ec- 
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Choose initial chemical potential and self-energy 



Calculate effective medium /^(w) 





Solve effective SIAM defined by ^, U 

and r(^-) 


Determine (n)siAM from solution 


Adjust /i to obtain desired (n) 



Iterate to self-consistency of S{z) 



Fig. 1. Flow diagram for the DMFT self-consistency cycle with 
fixed filling (n) . 



instead, with 

B(e,ui) = ——9m — = 

TT ^Uu,)C^{u,) - e 

and Cai'-^) = Lo+i-i — Sa-{^+iO^)- Obviously, expression (5) 
reduces to (4) without magnetic order, i.e. C<t('^) = Cai^^)- 
Throughout the rest of the paper we concentrate on 
results for a simple hypercubic lattice and 

f — t if ?-, j are nearest neighbors 

^ \ else ' 

The resulting DOS p^°\e) then becomes [10,12] 

p(°)(e) = -i=e-(^/**)^ . (6) 

t VTT 

In the following we use t* = = 1 as our unit of 

energy. 



tive medium, and the procedure is repeated until conver- 
gency is reached. It should be noted that, for a metastable 
state, no true convergency can be reached in the sense 
that the solution, when iterated further without adjusting 
the filling properly, will flow again into the phase sepa- 
rated one. Typically, for a converged calculation, the on- 
site energy between successive DMFT iterations shows a 
weakly damped oscillatory behavior but does not change 
by more than 2-3%. In order to minimize errors introduced 
by these oscillations we average all quantities over several 
iterations. 

To find the correct ground state, we need to calculate 
the ground state energy 



E 
N 
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U 



N 



{H} = -{H,} + -J2{n^^n,0 , (3) 



N 



where Ht is the kinetic part of the Hamiltonian (1). The 
expectation value {ni^ni^} can be determined within the 
NRG directly. The quantity (Ht), on the other hand, de- 
pends on the phase we are looking at. For the para- and 
ferromagnetic phases it is simply given by [11] 



^{Ht)=Y^ J rfeepW(e) J <L;f{uj)A^{e,u;) 



(4) 



with p^°^ (e) the density of states (DOS) for the non-inter- 
acting system, /(w) the Fermi function and 



1 



w + fi - € - + iO+) 



the spectral function of the Hubbard model in the DMFT, 
i.e. with k-independent one-particle self energy. 

For an antiferromagnetic state with Neel order one has 
to take into account the AB-lattice structure and the for- 
mula becomes [11] 

oo oo 

^{Ht)^2 J deepW(e) J dwf{u^)B{e,oj) (5) 



2.2 Weak-coupling results 

Let us briefly review some weak-coupling results as these 

will be frequently referred to in section 3. Since the hy- 
percubic lattice is a bipartite lattice, one obtains in low- 
est order, i.e. in Hartree approximation, a transition into 
a Neel state for any U > at T = below a critical 
doping 6^{U). For small J7 ^ the magnetization m as 
well as the critical doping depend non-analytically on U, 
i.e. m, 5^ oc exp (-l/(C//9(°)(0))) /[/ independent of the 
dimension. 

A quantity of particular interest in the DMFT is the 
single-particle Green function. The general structure of 
the Green function in the Neel phase for both Hartree 
theory and DMFT is given by expression (2) , where in the 
Hartree approximation Ea{z) reduces to S^{z) = Un^ = 
^U{n — am) with n the filling and m the magnetization. 
The local Green function is obtained from (2) by summing 
over k e MBZ, which yields for example for spin up 



g(°)(^C^HCm) (7) 



with Ccr(t^) = uj + iO"*" + IJ,— + a^m and 



oo 



de 



p(°)(e) 



(8) 



For the further discussion let us define 
= — /U — 



Then, as long as w < or w > a;+, the radicant in (7) is 
positive and the resulting DOS can be expressed as 
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For uji < u! < w+j on the other hand, the radicant in (7) 
is negative, i.e. -^C|(a;)Cx(a;) = i-\/|CT(^)Ci('^)l- Since for 
the particle-hole symmetric DOS (6) the Green function 
G^^\z) defined in (8) for purely imaginary arguments is 
purely imaginary, too, one finds 

i.e. the DOS has a gap between w- and a;+. As one ap- 



U|m| 



3 

<- 

Q. 




CO 
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Fig. 2. Behavior of the DOS in Hartree approximation close 
to the gap edges. 



to 

t5 




Fig. 3. The function $^{5/5i) from Eq. (11). Note the concave 
curvature between 5 = and 5 — \/25i . Tiie dashed line shows 
the actual behavior of the ground state energy following from 
a Maxwell construction. 

a" = 25?//9(°)(0). The function $^{5/5i) appearing in 

expression (10) leads to the full line in Fig. 3. Apparently, 
this function is not convex for small (5, i.e. the resulting 
phase is thermodynaniic;ally unstable towards phase sep- 
aration for dopings less than 5c = V^5i. The resulting 
ground state energy is then obtained from a Maxwell con- 
struction, given by the straight dashed line in Fig. 3. 



proaches oj- from below or a;+ from above, it is easy to 
confirm that 




Um 



(9) 



The corresponding DOS for a —I has a similar behavior. 
Here, however, the DOS diverges like — at the 

upper gap edge, and vanishes like i/|a; — a;_| at the lower 
one. 

In order to determine the thermodynamically stable 
phase one has to calculate the ground state energy as func- 
tion of the doping 6 = 1 — n. The result up to second order 
in U is [14] 



E{d) - E{0) = ~d + a" • ^"((5/(5i) 



where 



(If 1 

7:X 1 -X 

2 V 4 



I 4 



, 1 . 
1 + -X' 



X <1 



X>1 



(10) 



(11) 



and 6i is the critical doping for antiferromagnetism in 
Hartree approximation. The coefficient is given by 



3 Results 

Let us start with a short overview of the behaviour at 
half filling, n = 1. Here, the Neel phase is energetically 
stable. The variation of the DOS for increasing U from 
U — 1 (full curve) to [/ = 6 (dashed curve) is shown 
in Fig. 4. As expected, the DOS for small U resembles 
the form (9) predicted for weak-coupling, i.e. one sees the 
remnants of the characteristic square-root divergency in 
the spin up DOS at the lower gap edge and a correspond- 
ing power law at the upper gap edge. Tlicsc characteristic 
features however vanish rapidly with increasing U, and al- 
ready for = 3 the DOS mainly consists of the Hubbard 
peaks at Ll! = +U/2 and cj = —U/2 for a =1 and a =|, 
respectively; reminiscent of the behavior expected for the 
Mott-Hubbard insulator, where only the incoherent charge 
excitation peaks at high energies arc present [2,11,29,30]. 
Note that neither from the spectra in Fig. 4 nor from the 
behaviour of the magnetic moment in the inset of Fig. 4 
one can infer that at J7 w 4.1 the Mott-Hubbard metal- 
insulator transition occurs in the paramagnetic state [29, 
30]. 

Keeping U fixed at [/ = 3 and increasing 6 leads to the 
spectra shown in Fig. 5. Quite interestingly, the typical 

weak-coupling characteristics reappear in the spectra for 
small doping and are still recognizable for 6 = 13%. Note 
also that upon variation of doping and hence of the mag- 
netization the spectra are not shifted in the same way as 
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Fig. 4. DOS for spin up and down at half filling in the anti- 
ferromagnetic phase as function of U . While for small values 
of U the weak-coupling form (9) is approximately reproduced, 
the DOS for large U is basically that of the Mott-Hubbard 
insulator. The inset shows the magnetization as function of U . 



in Hartree theory. Instead, the dominant effect is a strong 
redistribution of spectral weight from the Hubbard bands 
to the Fermi level. Eventually, in the paramagnetic phase 
one recovers the well-known three peak structure of the 
doped Hubbard model in the DMFT [11]. 

The evolution of the spectra both at and off half filling 
can be understood within a simple picture. In Fig. 6 we 
show a sketch of the Hartree bandstructure of the Hub- 
bard model in the Neel state, which has two branches in 
the MBZ and a gap of width oc U\m\ between them. If, 
on the other hand, we inspect the paramagnetic solution, 
one for example finds at half filling and for small values 
of f/ a Fermi liquid with quasiparticles defined on an en- 
ergy scale larger than f7|m|. This situation is indicated 
by the arrow labeled (a) on the left side of Fig. 6. Here 
we expect, and indeed find for the antiferromagnetic so- 
lution (see full curve in Fig. 4), a DOS that shows the 
characteristic van-Hove singularities of Fig. 2. Increasing 
U eventually leads to a situation, where the energy scale 
for the quasiparticles in the paramagnetic state is finite 
but much smaller than U\m\ (arrow (b) in Fig 6). The 
self-energy in the energy region of the van-Hove singular- 
ities then has a large imaginary part and will completely 
smear out the characteristic structures. Further increasing 
U into the Mott-Hubbard insulator will then not change 
the picture qualitatively, explaining the similarity between 
the curves for ?7 = 3 < Umyt and U = &> C/mit in Fig- 4. 
With finite doping, we move the chemical potential into, 
e.g., the lower band; this means that even for a relatively 
small quasiparticle energy scale one again sees the van- 
Hove singularities at the band edge, which results in the 



Fig. 5. DOS for spin up and down for [/ = 3 and different 
dopings 5 = 7%, 5 = 13%, 5 = 16% and 5 = 20%. The system 
at 5 — 20% is already in the paramagnetic phase. 







(a) 
















1 





cU|m| 



-1/2 





k/7t 



1/2 



Fig. 6. Schematic picture of the Hartree bandstructure of the 
Hubbard model in the Neel state. The arrows at the left hand 
side of the figure represent the energy scales of the correspond- 
ing paramagnetic Fermi liquid for half filling and weak coupling 
(a), half filling and intermediate coupling (b) and finite doping 
and intermediate coupling (c). 



well defined singularities in the spectra for small doping 
in Fig. 5. 

From the occupation numbers obtained after con- 
vergence of the DMFT calculation one can calculate the 
magnetization per electron, m — [n^ — n|)/(ri| -I- n^), as 
function of the doping 5. The results for U — \ and U — ?> 
are shown in Fig. 7a and b together with a fit to a power 
law 

5 



m{S) — mo 



1 - 



(12) 
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Fig. 7. Doping dependence of the magnetization per electron 
for (7 = 1 (a) and U = 3 (b). The full lines are fits with the 
function (12), the resulting fit parameters are summarized in 
Table 1. 



The resulting fit parameters are summarized in Table 1. 
As expected for a mean-field theory, the value for the crit- 
ical exponent is v = 1/2. 




Fig. 8. Doping dependence of the energy of the paramagnetic 
phase (squares) and the Neel phase (circles) for U — 1 (a) and 
U — 3 (b). The full lines are fits with the function (13), the 
corresponding fit parameters are summarized in Table 1. The 
dashed lines are the result of a Maxwell construction for the 
ground state energy. 



Finally, with the converg ed DMFT self-energy E^{z) 
we can calculate the expectation value {H)/N according 
to Eqs. (3) and (4) respectively (5) for the paramagnetic 
and antiferromagnetic phase. The results for the charac- 



teristic function 

<?(<5)-i?(<5) + |5-i?.„ag(0) 

are summarized in Fig. 8a and b. In Fig. 8 the energies of 
the antiferromagnetic phase are represented by the circles, 
those of the paramagnetic phase by squares. The full lines 
interpolating the antiferromagnetic data are fits to the 
function 

<Z>(<5) = a<Z>H(V<5i) + 7(^) (13) 

with <?^(x) according to (11). The fit parameters are sum- 
marized in Table 1. The use of the function <?^(x) in (13) 



Magnetization Energy 
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mo 




u 




&i 
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1 


0.4 


0.06 


0.49 


0.07 


0.047 


0.52 





3 


0.9 


0.16 


0.54 


0.157 


0.191 


0.33 


0.026 



Table 1. Results of the fits of m(5) in Fig. 7 to expression 
(12) and E(S) in Fig. 8 to (13). 



is motivated by the results of van Dongcn [14]. The lines 
interpolating the paramagnetic data are meant as guides 
to the eye only. The dotted vertical lines denote the value 
(5^^ as obtained from Fig. 7. 

The antiferromagnet obviously has the lower energy 
as compared to the paramagnet in the region < (5 < 
(5^^. However, in both cases [/ = 1 and J7 = 3 we find a 
clear non-convex behavior in E(S) in that region, i.e. the 
aforementioned signature of an instability towards phase 
separation. The true ground state energy as function of 
b is obtained again via a Maxwell construction, leading 
to the dashed lines in Fig. 8 and the values i5^^ given in 
Table 1. Note that in both cases « (5^® within the 
accuracy of the fitting procedure. 

While for U ~ 1 the function 'P{b) nicely follows the 
weak-coupling prediction (10) with renormalized constant 
a one finds a sizeable contribution ~ (5^ for J7 = 3. This 
additional term results in a much weaker non-convex be- 
havior of E{b) for [/ = 3. 

For values [/ > 4 we were not able to find a stable 
solution with Neel order and well-defined doping S > 0, 
although for b < bc{U) the paramagnetic phase becomes 
unstable. However, the numerical calculations rather pro- 
duce a cycle encompassing a range of fillings instead of 
one solution with definite filling here. It might be inter- 
esting to note that at least each of the fillings in this cycle 
has a unique magnetization associated with it and that all 
spectra in this cycle correspond to an insulator. Currently 
it is neither clear what type of magnetic solution we find 
here, nor whether the breakdown of the Neel solution is 
a true physical effect or due to numerical problems. Since 
at half filling the Neel state is present at these values of 
U, incommensurate structures or again a phase separated 
state seem to be possible. 
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For values of U beyond Uc ~ 25 yet another magnetic 
phase appears, namely the ferromagnet. The existence of 
this phase has been observed in the case of a hypercubic 
lattice and U = oc' [31] and for a generalized fee lattice 
[32] before. Since these calculations had to be done at 
finite and comparatively high temperatures, questions re- 
garding the ground state magnetization and, especially in 
the case of a hypercubic lattice, the actual extent of the 
ferromagnetic phase in {S, U) space could not be discussed 
satisfactorily. 

As an example for the ferromagnetic phase at T = 
Fig. 9 shows the magnetization per electron, m{5) = (n| — 
n|)/(n| +ni) as function of doping (Fig. 9a) and the local 
DOS for two dopings (Fig. 9b) {U = 50). The data for 
m{5) in Fig. 9a are fitted to the function 



m{6) = Too • 




(14) 



and the result is given by the dotted line. The parame- 
ters for the fit are too = 1, (5c = 14.6% and v = 2.75. 
While for 5 y Sc the typical mean-field behavior, i.e. 

m{S) oc Y^l — ^, is obtained, the result for 5 — > is 

/ s 2.75 

rather unconventional, viz m{6) oc 1 — i I j- ) . This 



fit assumes that a fully polarized state is only reached as 
6^0 [33]. Note, however, that the numerical results for 
the magnetization m{S) for small 5 are also consistent with 
a fully polarized ferromagnet at finite S. 

It is also quite apparent from the DOS in Fig. 9b, that 
the ferromagnetism found here cannot be understood on 
the basis of the typical Stoner theory. In contrast to the 
shifts of the spectrum expected in the latter, we observe a 
strong redistribution of spectral weight instead, but retain 
otherwise the typical structures due to the strong correla- 
tions. Only in the case 5—^0 the spectrum again resembles 
that of a free system for the (almost completely polar- 
ized) majority spins. The minority spins become strongly 
depleted below the Fermi energy, the spectral weight can 
be found almost completely in the upper Hubbard band 
situated around lo Ki U/2 (not shown in the figure). Never- 
theless, we observe a tiny resonance just above the Fermi 
energy even as 5 ^ 0. 




0,5 




0,5 



— ^ 





1 (b) 


1 1 1 1 

— 5=12% 

- - -- 6=1.7% 


, ■ | — 



-3 



-2 



-1 



CO 



Fig. 9. (a) Ferromagnetic magnetization per electron as func- 
tion of doping 5 ioT U = 50. The full line is a fit to the function 
(14). The critical doping is Sa « 14.6%. Note that for 5 
the results arc consistent with a fully polarized ferromagnetic 
state, (b) Local density of states for U — 50 and two character- 
istic dopings 5 = 12% (full linos) and S = 1.7% (dashed linos). 
In contrast to the Stoner theory, one finds comparatively small 
shifts in the spectra, but a strong redistribution of spectral 
weight. 



4 Summary and conclusion 

In this paper, we used the dynamical mean-field theory 
together with Wilson's numerical renormalization group 
to investigate the ground-state properties of the Hubbard 
model on a hypercubic lattice with nearest-neighbor hop- 
ping both at and off half-filling. While at half-filling the 
ground-state is antiferromagnetic for all t/ > 0, at least for 
the weak and intermediate coupling regime this magnetic 
order can only be realized in a phase-separated state for 
any finite doping, thus supporting and extending earlier 
weak-coupling predictions. 

The mapping of the Hubbard model for large U to an 
antiferromagnetic t-J model strongly suggests the domi- 



nance of antiferromagnetism in the ground state. The re- 
sults for the Hubbard model in this paper show, however, 

that the type of magnetic order for intermediate values of 
the Coulomb repulsion U off half filling is still an open 
issue; furthermore, the role of phase-separation (which is 
observed for U < 3) has still to be clarified for larger val- 
ues of U. 

The results are summarized in the schematic {S, U) 
ground state phase diagram of Fig. 10. To allow the inclu- 
sion of all values < U < oo, the ratio U/ (1 -|- U) is used 
on the abcissa. Close to half filling we find a phase sepa- 
rated Neel antiferromagnet up to a certain value oiU < 4. 
The magnetization as function of doping follows a typical 
mean-field behavior in all cases studied and the spectra 
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Fig. 10. Scliematic ground state phase diagram. At half-filling 
{5 — 0), the ground state is antiferromagnetic. Close to half fill- 
ing and small U, a phase separated antiferromagnet is realized. 
For values of U beyond U ~ 4 a. magnetic phase is observed, 
whose precise character could however not be identified. For 
very large U > 25 and up to 5 « 30% a ferromagnet is found. 

show the characteristic van-Hove singularities of the band 
structure in the Neel state in cases where the characteris- 
tic energy scale of the paramagnet is large enough. Most 
important is the observation that, as typical for correlated 
ordered systems, the spectra are not strongly shifted, as 
e.g. predicted by Hartree theory, but rather show a strong 
redistribution of spectral weight. 

For values of /7 > 4 the system shows the tendency 
towards a magnetic instability, which could not be further 
identified due to technical problems in the solution of the 
DMFT self-consistency. However, we at least can exclude 
ferromagnetism here and a speculative possibility would 
be the occurence of incommensurate phases or magnetic 
phases with additional charge order. While the former can- 
not be addressed easily within the present method, the 
latter possibility will be investigated further. 

At very large values of t7 > 25, there is a region of fer- 
romagnetism, extending between < S < 30% as J7 ^ oo 
[31] . For a fixed value of U, the magnetization per electron 
in the ferromagnetic state shows the tendency to saturate; 
from the numerical data it is of course impossible to re- 
liably conclude whether the system is fully polarized at 
a finite 6 already or only as (5 ^ 0. The data are consis- 
tent with both scenarios, but the latter is supported by 
analytical treatments of the case (5 — > 0. As in the case of 
the antiferromagnet, the spectrum shows a rather strong 
redistribution of spectral weight, not simply a shift of the 
features, as would be expected from Stoner theory. 

The phase diagram shows a peculiarity, which has al- 
read been pointed out by Obermeier et al. [31]. In the 
region of very large U and (5 — > there exists the possibil- 
ity of a direct transition between the "antiferromagnetic" 
phase and the ferromagnet. As at the point {S, U) = (0, oo) 
all possible spin configurations are degenerate, one can 
speculate how the phase diagram looks like as {5, U) ~* 
(0,oo) [34]. Generic possibilities are sketched in Fig. 11. 
There can either be a direct transition between the two 



(O,--) ^5 (0,~) ^5 (0,~) 




Fig. 11. Possible realizations of the phase diagram as {5, U) 
(0, oo): A direct transition between an antiferromagnet and a 
ferromagnet as in (a) , a small paramagnetic phase between the 
two las in (b) or a mixed type of phase (e.g. ferrimagnet) as in 
(c). 

phases (Fig. 11a), which quite likely would then be of first 
order, a gap filled by a paramagnetic phase (Fig. lib) or 
a new phase, e.g. a ferrimagnet interpolating between the 
two extremes. The a priori exclusion or verification of any 
of these structures is, without a detailed knowledge of the 
analytic behaviour of the relevant quantities as function 
of ((5, U) in the vicinity of (5, U) = (0, oo), not possible. 

While there is a consensus about the magnetic prop- 
erties of the Hubbard model in a qualitative sense, the 
direct inspection of details still reveals unexpected sur- 
prises. Even within the DMFT, where one can safely state 
that the paramagnetic phase diagram including the Mott- 
Hubbard metal insulator transition is now understood, the 
investigation of the magnetic properties is far from com- 
plete. Obvious open questions are the magnetic properties 
at intermediate values of U and the behavior when the an- 
tiferromagnetic and ferromagnetic phases meet. Further- 
more, the behaviour of the magnetic phases, especially the 
spectral properties in the ordered phases, in the presence 
of a frustration due to longer range hopping, has not been 
addressed yet. This might be of some interest regarding 
the question how the first order Mott-Hubbard transition 
manifests itself in the magnetically ordered state. Work 
along these lines is in progress. 

We acknowledge useful discussions with P.G.J, van Dongen, 
M. JarreU, J. Keller, D. Logan, M. Vojta and D. VoUhardt. 
This work was in part supported by the SFB 484 Kooperative 
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